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We identify a new kind of elementary excitations, spin-rotons, in tfie doped Mott insulator. They 
play a central role in deciding the superconducting transition temperature Tc, resulting in a simple 
Tc formula, fcaTc — Eg/G, with Eg as the characteristic energy scale of the spin rotons. We show 
that the degenerate S = 1 and S — rotons can be probed by neutron scattering and Raman 
scattering measurements, respectively, in good agreement with the magnetic resonancelike mode 
and the Raman Aig mode observed in the high- Tc cuprates. 

PACS numbers: 74.20.Mn, 71.10.Hf, 71.10.Li 



I. INTRODUCTION 

To fully understand the nature of high- Tc supercon- 
ductivity in the cuprates, one essential task is to identify 
the most important elementary excitation which controls 
the superconducting transition. 

In a conventional BCS superconductor, the Bogoliubov 
quasiparticle constitutes the most crucial low-lying exci- 
tation. In a d-wave state, nodal quasiparticle excitations 
generally lead to a linear-temperature reduction of the 
superfluid stiffness Ps 

Ps{T) = ps{Q) - aT (1) 

which, however, would be normally extrapolated to a 
transition temperature {ps(Tc) = 0) much higher than 
the factual Tc in the cuprates, based on the microwave 
measurements of the penetration depth which determines 
the superfluid density^. 

On the other hand, in view of the small superfluid den- 
sity in the cuprates, which are widely considered to be a 
doped Mott insulator—, the phase fluctuation of the su- 
perconducting order parameter has been suggested^ to 
play an important role in the transition regime, which 
can be characterized by the following London action 

L^^ Jd'ri^^ + qA^r (2) 

where specifies the U{1) phase of the order parameter 
of condensate carrying charge q, and A"^ is the external 
electromagnetic field. In this point of view, the supercon- 
ducting transition is of a Kosterlitz-Thouless (KT) typo^ 
with the proliferation of topological vortices 

^ dr • V(/) = ±27r (3) 

which destroy the phase coherence of superconductivity 
resulting in ksTc ~ PsiT~). 

However, a striking and puzzling empirical Tc for- 
mula for the cuprate superconductors has been known 
experimentally^i^ '^^i^^ , which is simply given by 

ksTc = ^ (4) 
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FIG. 1: [Color online] The characteristic energies observed by 
inelastic neutron scattering (INS) and electron Raman scat- 
tering (ERS) experiments versus the superconducting transi- 
tion temperature Tc for the high- Tc cuprates. The straight 
line shows the empirical formula ([4]), which will be derived 
in the present work. Here the solid squares represent the 
INS resonance mode, with different colors indicating different 
families including hole-doped ¥123^^, Bi2213i4, T1220li^ and 
IIgl201— , and electron doped Pro.88LaCeo.i2Cu04_,r— ; the 
solid circles represent the ERS Aig mode, including the hole 
doped Y123i^i^, Bi2213^, T12201— i^^, Hgl201-^, La2142^, 
T122132^, T122232^, Hgl2132£, Hgl2232i and the electron 
doped NCCO compound^^. 



where k ^ 6 and Eg denotes the charcteristic en- 
ergy scales observed in inelastic neutron scattering 
(INS)li^i^ii^iHiiia5ii6ii2 and electronic Raman scattering 
(ERS)ia'i^^'20^'22^'24^^^^ measurements, as il- 
lustrated in Fig. [TJ Here Eg in INS corresponds to the 
well-known resonance energy^ in the literature, which is 
a spin-triplet excitation at momentum centered around 
the AF wave vector Qaf = (T^jTr). By contrast, Eg in 
ERS corresponds to a singlet mode in the Aig channel 
near momentum Qo = (0,0). The ERS data in Big and 
B2g channels have provided the compelling evidence for 
the d-wave pairing symmetry in the cuprate superconduc- 
tors, however, the Aig peak at Eg remains an unresolved 
mysterjii^. As shown in Fig. [1] more materials can be ac- 
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cessiblc by ERS than INS, including the La2_a;Sr2;Cu04 
compound in which there is no direct INS evidence for a 
sharp resonancehke mode but a singlet mode in ERS^ 
has been still found with Eg well fit by Q. 

The above empirical scaling law of vs. Eg implies 
that the elementary excitations controlling the supercon- 
ducting transition in the cuprates should be composed of 
two degenerate modes, with quantum number S = and 
1, respectively, as probed in ERS and INS. Note that in 
the literature the magnetic resonancelike mode observed 
in INS has been sometimes interpreted as the bound state 
of a Bogoliubov quasiparticle pair near the antinodal 
regime due to the residual superexchange interactiorti. 
In this picture it would be hard to understand the neces- 
sity for the existence of a singlet bound state with the 
roughly degenerate energy. The further challenging and 
fundamental question is, given the presence of two de- 
generate modes observed in ERS and INS, how can they 
directly influence the superconducting coherence? 

A proposal made by UemuraF^ recently is that the two 
quasi-degenerate modes observed in INS and ERS may 
originate from soft modes in spin and charge channels 
in an incommensurate stripe state, which are called^i 
twin spin/charge roton mode, in analogy with the soft 
phonon-roton mode towards solidification in superfluid 
*He. Hence the mechanism for superconducting transi- 
tion is due to the substantial reduction of the superfluid 
density by thermal excitations of such twin spin-charge 
soft mode at Tc/2 < T < Tc, whereas the quasipar- 
ticle excitations mainly dominate at lower temperature 
< Tc/2. 

Nevertheless, according to the experimental results 
shown in Fig. [1] it seems that the Tc formula (|4]) holds 
more generally than simply in a neighborhood of stripe 
states^^. It calls for an intrinsic "spin-charge entan- 
glement" in the superconducting phase of the cuprates. 
Namely, magnetic excitations at Qaf should have some 
kind of profound effect on the superconducting conden- 
sation such that thermal excitations of the former can 
be destructive to the latter, much more effective than 
the usual nodal quasiparticles in the BCS theory. Fur- 
thermore, the mechanism should allow for a degenerate 
singlet mode, which may be not associated with a soft 
mode of any charge order as its characteristic momen- 
tum is around Qo, to play an equally important role. 
Lastly, the simple scaling relation ^ with a universal 
K should be independent of the details of materials in- 
cluding the charge inhomogeneity. Or more precisely, all 
the detailed properties of the system should influence Tc 
mainly through the characteristic energy scale Eg. 

In this paper, we will demonstrate that a self-consistent 
mathematical description of superconductivity in doped 
Mott insulators can give rise to a systematic account for 
the above-mentioned novel properties including the Tc 
formula ([4]). In the superconducting state, besides the 
emergent quasiparticles as the recombination of charge 
and spin, the most nontrivial elementary excitations are 
the vortex-antivortcx bound pairs locking with free spins 




FIG. 2: [Color online] Schematic illustration of an S = 
(singlet) and an S = 1 (triplet) spin-rotons. Each of them 
is composed of a supercurrent vortex-antivortex bound pair, 
with a pair of neutral free spins sitting at the two poles of 
the two-dimensional roton. Such a spin-roton composite is an 
elementary excitation in the superconducting state of a doped 
Mott insulator described by the phase string theory^^^. 



at the poles, with total spin 5 = or 1, as illustrated in 
Fig. [21 We shall call these excitations spin-rotons in the 
following, which are distinguished from those proposed 
by Uemurapii as they are not slaved with any charge and 
spin orders, but a direct consequence of the phase string 
effect^i^ in the t-J model with a peculiar nonlocal spin- 
charge entanglement: neutral spins locking with charge 
supcrcurrents22iS. 

These spin-rotons will naturally include two degener- 
ate excitations. The degeneracy of these modes with 
spin quantum number S = and 1 is due to the fact 
that the pair of neutral spins are excited "spinous" from 
an underlying rcsonating-valence-bond (RVB) spin back- 
ground. The degenerate spin-roton modes thus indicate 
spin-charge separation, but with a twist. That is, a sta- 
ble spin-roton object in the superconducting phase also 
implies a spinon- confinement as two spinous cannot be 
separated freely in space due to the logarithmic poten- 
tial between the vortex and antivortex. Such rotonlike 
supercurrents will play a central role in deciding the su- 
perconducting phase coherence transition as in ([1]). We 
will show that the singlet and triplet spin-rotons can be 
indeed directly probed by ERS in Aig channel at Qo 
and INS near Qaf- They have the minimal character- 
istic energy Eg ~ SJ in the low-doping regime with the 
magnitude in good agreement with the experimental data 
where S denotes the doping concentration and J is the 
superexchange coupling constant determined in the un- 
doped case. 

The remainder of the paper is arranged as follows. In 
Sec. Ill Al we introduce the description of a doped-Mott- 
insulator superconductor, obtained prcviously^^iS based 
on the phase string theory^Si^ of the t-J model, by using 
a phenomenological construction. We argue that in order 
to incorporate the influence of spin degrees of freedom 
(which is important in a lightly doped Mott insulator 
where spins constitute the majority of low-lying degrees 
of freedom) under the requirement of no time-reversal 
and spin rotational symmetry breakings, one is naturally 
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led to a modified action for superconductivity. In Sec. 
IIIBl tlie spin-roton excitations as a direct consequence of 
this formulation are discussed. Then, how the spin-roton 
excitations as the resonancelike modes can be probed by 
INS and ERS are discussed. In Sec. IlICi the Tc formula 
(III) determined by the spin-roton excitations is obtained. 
Finally, in Sec. IIIIl a general discussion will be given. 

II. SPIN-ROTON EXCITATIONS 

A. Phenomenological description of a 
doped- Mott-insulator superconductor 

Prom a doped Mott insulator point of view^, the su- 
perconductivity in the cuprates occurs in a small doping 
regime where the charge carrier number is greatly re- 
duced as compared to the total electron number. Namely, 
the strong on-site Coulomb interaction will make the 
charge degrees of freedom partially frozen, while the full 
spin degrees of freedom of the electrons remain at low 
energy. Thus the London action ([2]) should be modified 
in order to properly reflect the Mott physics. 

For example, in the U(l) slave-boson gauge theory 
description^, the charge carriers are described by spin- 
less bosons known as holons. The superconducting state 
corresponds to the Bose condensation of the holons, with 
([2]) replaced by 

i,. = ^ y d'r(V(/) + + eA^y (5) 

where ps is proportional to the density of condensed 
holons and q = +e, in contrast to the conventional Lon- 
don action where the condensate of Cooper pairs of the 
electrons is involved with q = — 2e. As a component of 
the electron fractionalization, holons are no longer gauge 
neutral and are generally coupled to an internal emergent 
gauge field A''. In the U(l) slave-boson gauge theory^i^ 
A** will be also minimally coupled to the other component 
of the electron fractionalization, i.e., neutral spins called 
spinous. However, since the latter are in RVB pairing, 
the internal gauge field A^ is expected to be suppressed 
due to the "Meissner effect" of the RVB state, whose 
mean-field transition temperature is presumably much 
higher at low doping. Consequently in such a mean-field 
"pseudogap" regime A** gains a mass and cannot play a 
role as a new source to effectively reduce T^. 

However, the U(l) slave-boson gauge theory is not 
the only possible theoretical description for the doped 
Mott insulator. In the following, we shall elucidate in a 
phenomenological way an alternative self-consistent con- 
struction. It will reveal the existence of a new math- 
ematical structure^i^^, in which the charge condensate 
can become strongly correlated with spin excitations. 

The key distinction will be that, instead of minimally 
coupling to both the holon and spinon currents in the 
U(l) slave-boson gauge theory, here A" will only mini- 
mally couple to the holon matter field as given in ^ , not 



to spinon currents. Instead its strength will be generated 
from the spinon matter field according to the following 
gauge-invariant relation 

(fdr-A' (r) = cj)o [ d?v [n\ (r) - n\ (r)] (6) 

Jc JSc 

Here the flux of A" within an arbitrary loop c on the left- 
hand-side (l.h.s.) is contributed by ±00 flux-tubes bound 
to individual spinous on the right-hand-side (r.h.s.), with 
n||(r) denoting the local density of spinous where the 
integration runs over the area Ec enclosed by c. 

Due to the sign change between the ] and [ spins on 
the r.h.s. of A''(r) will explicitly preserve the time- 
reversal (TR) symmetry, as t <-> | under the TR trans- 
formation. This is in contrast to a conventional elec- 
tromagnetic vector potential A'^, which breaks the TR 
symmetry. However, since the path c is oriented, the 
spin rotational symmetry may be broken for a general 
<j)Q. But under a specific choice 

fZ-O = TT (7) 

one finds that the spin rotational symmetry can be still 
maintained: without loss of generality, one can con- 
sider a loop c which encloses a single spin such that 
§^dY ■ A* = ±00 = ±7r which is still spin-dependent. 
However, such a spin-dependence sign change can be ef- 
fectively compensated in ([5]) by combining with a proper 
topological vortex of the holon condensate given in ([3]). 
Such a "large" gauge transformation will not cost any 
energy in ^ when A'^ = 0. It is also "legal" to precisely 
bind such a holon vortex core of ^ with the spinon be- 
cause the no double occupancy constraint in the doped 
Mott insulator dictates that a site without a holon must 
be alway occupied by a neutral spin. 

Hence, based on some general physical considerations, 
the London action for a superconducting state can be 
modified in a fundamental way in a doped Mott insula- 
tor, with an internal vector potential A" emerging as a 
topological gauge field without breaking the time-reversal 
and spin rotational symmetries. 

According to (O, the charge current will be determined 
by the London equation 

J,, =p,(V(/) + A^ + eA^) (8) 

For an isolated neutral spin, in terms of ([6]), there will be 
vortexlike charge currents induced from the charge con- 
densate with § dr-Jh = ±Ps7i' in the absence of A*^, where 
± will be independent of the spin index based on the 
above discussion. Namely each neutral spin can induce a 
current vortex with two opposite vorticities as illustrated 
in Fig. [31 which is known as a spinon-vortes^2i2^. 

According to a general argument given by Haldane and 
Wu^, since a spinon behaves like a supercurrent vor- 
tex, its motion through a closed path c must then pick 
up a Berry's phase which is determined by the number 
of superfluid particles of the condensate in the area Sc 
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FIG. 3: [Color online] Schematic illustration of single spinon- 
vortices. An isolated neutral spin (spinon) in the supercon- 
ducting state will always induce a vortexlike supercurrent re- 
sponse from the charge condensate according to the gener- 
alized London action ([5}. Notice that the vorticity sign of 
the vortex is actually independent of the spin orientation as 
long as (/!>o = TT in (O, which preserves the spin rotational 
symmetry. 

enclosed by c, as if it sees an effective "magnetic-field" 
described by a vector potential A'': 

A$Borry(c) = 'Po J /Oft (r) 

= <j)Jr-A^{r) (9) 

Here p/i(r) denotes the local supcrfluid density of con- 
densed holons, with (j)Q = n instead of 2tt. 

Thus, one may write down a minimal gauge- invariant 
Hamiltonian for spinous simply as 

= -Js C^J-.e'"^*^ (10) 

<ij>(y 

where defines the bosonic creation operator for a 
spinon at site i with a spin index a. Here is the 
lattice version of the gauge potential A'' (r) introduced 
in ^ and the sign a in front of the gauge phase in (fTUl) 
will ensure the TR invariancc. 

Although one can alternatively write down an effective 
model with the hopping term bl^bjae'''^^'^ replacing the 
RVB pairing term 

C^j-aS*'^'^''' ill Cni)! without break- 
ing the gauge and TR symmetries, (jlOp is physically more 
meaningful because in the ground state spinous will be 

all paired up with (bl^b]_^e''"^'^^ ^ = A" /2 ^ which au- 
tomatically satisfies the spinon-confinement requirement 
to ensure superconducting phase coherence as to be dis- 
cussed below. Furthermore, dr ■ A'' = at half-filling, 
where Hs PH)) reduces to the Schwingcr-boson mean-field 
Hamiltonian which well captures the antiferromagnetic 
(AF) correlations including the long-range AF order at 
T = 0^. 

Therefore, the London action ^ for superconductiv- 
ity has been phenomenologically modified for the doped 
Mott insulator in ([5]). Here the charge condensate will 



be generally coupled to neutral spin excitations, ubiqui- 
tously presented in a doped Mott insulator governed by 
PU)) . via an emergent topological gauge field Such a 
self-consistent description based on ([5]), (O, Q and pU]) 
can be justifie d? -- -" in the phase string theory of the t- 
J model, with the superfiuid stiffness = ph/rrih {m^ 
is the effective mass for holons) and effective coupling 
constant Js in (1 101) determined microscopically. One is 
referred to Ref. [SJ and the references therein for details. 
Although it is not a unique construction for a doped Mott 
insulator (one can alternatively have other possible math- 
ematical constructions like the U(l) slave-boson gauge 
theory description^, for example, as mentioned before), 
some very unique consequences will follow from such a 
self-consistent approach, which can be directly compared 
with experiments. 



B. Spin-roton excitations 

A direct physical consequence is that a single spinon 
excitation in the superconducting state will not be per- 
mitted because the self energy of a vortex shown in Fig. 
[3]is logarithmically divergent. Then all the spinous in the 
superconducting state must appear in pairs, with the as- 
sociated supercurrent vortices forming vortex-antivortex 
bound pairs, as illustrated in Fig. [51 These bound ob- 
jects are referred to as spin-rotons, which carry total spin 
(singlet) and 1 (triplet), charge 0, together with a su- 
percurrent structure analogous to a two-dimensional ro- 
ton excitation in a Bosc condensate. In this sense, the 
spinous must be "confined" and only integer spin excita- 
tions are allowed in the superconducting state. 



1. Resonancelike characteristic energy Eg 

The spinon Hamiltonian (jlOp can be easily 
diagonalizcd^l under the condition that the holons 
are uniformly condensed with ph = Sa~^ {a is the lattice 
constant) as outlined in Appendix [X] The solution 
of PU)) has an uneven Landau-level-like spectrum for 
spinon excitations as shown in the inset of Fig. [4l which 
arc excited by breaking up RVB pairs in the ground 
state. 

At low temperature, we shall focus on the lowest ex- 
cited level at Es = Eg/2 for simplicity. In the main panel 
of Fig. \^ Eg = 2Es is shown as a function of doping un- 
der a proper consideration of the doping dependence of 
J^. The corresponding spinon wavepacket looks like 

with a "cyclotron length" Oc = a/Vrr^. Namely, the 
lowest spinon excitations governed by pop are non- 
propagating modes of an intrinsic size in order of Oc. 
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FIG. 4: The doping dependence of the characteristic energy 
scale Eg of the spin rotons is shown. Here Eg = 2Es with Es 
as the lowest excited energy of the spinon spectrum shown in 
the inset, obtained based on pO|) at a specific doping concen- 
tration 6 = 0.125. 

Here the degenerate levels are labeled by the coordi- 
nates Rmr^, the eentcrs of the spinon wavepacket (fTTjl . 
which form a von Neumann lattice with a lattice constant 
^0 = VSTrCc, as shown in Fig. [5] 

After integration over the original lattice index in 
the modified London action ([5]) at A"^ = 0, one can obtain 
(see Appendix[B| an effective interaction term for spinon- 
vortices on the von Neumann lattice 

TT 1 |I^m^R-m'| 

t'int = >^ In q-mq-m' (12) 

4 — CD 

where q^, (= ±1 or 0) denotes the vorticity for each 
spinon- vortex on the site Rm, and to avoid the logarith- 
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FIG. 5: [Color online] The degenerate spinon modes in the 
lowest energy level, shown in the inset of Fig. [l] are labeld 
by Rm which form the von Neumann lattice with a lattice 
constant = V^clc with the cyclotron length Oc — u/Vt^ 
as the size of each spinon wavepacket. Here the case 5 = 1/8 
and ^0 = 4a is shown. For each Rm, there is an additional 
degeneracy 3 = 4, corresponding to orthogonal wavefunctions: 
Wmtiri), Wmi(rO, {-iy'Wm^-{ri) and {-ly^w^iiri)^. 



mical divergence, the neutral constraint ^^^^ q„i = will 
be imposed. So the total energy of the spinon-vortices is 
given by 

m 

It is noted that there is a four-fold degeneracy, g = 4, at 
each site Rm as mentioned in the caption of Fig. [5] 

Note that a conventional vortex-antivortex pair in a 
KT system will normally shrink at low temperature and 
be annihilated in the ground state. But a spin-roton in 
the present case cannot literally disappear in the ground 
state because the two spins sitting at the poles of a roton 
in Fig. [2] cannot annihilate each other. Nevertheless, the 
roton supercurrents surrounding the neutral spins will 
have minimal effect on the ground state. In fact, as the 
solution of pO|) . spins will form short-range RVB pairs 
in the ground state, of a characteristic length scale ~ ac 
which is comparable to the finite core size of each pole 
of a spin-roton in Fig. [2] (the spin trapped at the core 
cannot sit still due to the uncertainty principle and the 
intrinsic core size is set by the cyclotron length Oc). Thus, 
the surrounding rotonlike supercurrents around an RVB 
pair will be effectively canceled out in the ground state. 
In other words, the London action ^ will be decoupled 
from the neutral RVB spin background as A" w and the 
excited spinon-vortices are effectively described by (jl3p . 

Hence the spin-roton structure shown in Fig. [2] will 
emerge as the bound pair of the excited spinous, which 
are of spin triplet {S = 1) and singlet {S = 0), respec- 
tively, and degenerate in energy. The spin-rotons here 
will have a minimal energy scale Eg when two excited 
spinous are located at the same von Neumann lattice site 
such that the vortex-antivortex supercurrent structure is 
effectively annihilated with Uint = 0. 

The degenerate singlet and triplet spin-rotons imply 
the spin-charge separation: i.e., the existence of single 
spinous carrying = 1/2 and zero charge as individual 
excitations, which do not interact with each other mag- 
netically. However, we have also seen that these spinons 
must be confined spatially in pairs, appearing at the poles 
of roton supercurrent structure and subjected to loga- 
rithmic attraction Uint- Therefore, in such a non-BCS 
superconducting state the spin-charge separation has a 
twist, which is characterized by new elementary exci- 
tations of degenerate spin-rotons instead of individual 
spinons. In other words, the spinon confinement does not 
mean a spin-charge tight recombination like in a conven- 
tional Fermi liquid or BCS superconductor of the elec- 
trons. Rather, at a short distance scale ~ ^o- the confin- 
ing force Uint becomes negligible and the spinons are still 
"asymptotically free" . 

2. INS and EES probes 

Experimentally, the neutron and Raman scattering 
measurements can provide direct means to probe such 
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novel excitations, in spin triplet and singlet channels, re- 
spectively. 

Define the spin-spin correlation function 

X..(T,r,-r,) = -(TA^(r)5|(0)) (14) 

where r denotes the imaginary time, 5f = ^ '^^\a^i'^- 
Similarly a density-density correlation function which 
can be detected by the electron Raman scattering^ is 
defined as follows 

XERS = -(r.rA„(r)TA,,(0)) (15) 

where the Aig density operator— tai^ = 

-5 S{ij>a cLcjo- + h.c. Here is the electron op- 
erator whose relation with the holon and spinon 
operators is given in Appendix [XI 

Based on the Bogolivbov transformation (jAl|) and the 
phase string representation for the electron operatore da , 
one can express and Ci^ in terms of and ^ma as 
shown in Appendix A. We shall mainly concentrate on 
energies near the minimal Eg^ where the total Hamilto- 
nian reduces to (|13p in which the interaction term [/jnt 
can be also neglected because Sf and tai^ only create a 
pair of spinons locally within a von Neumann lattice site 
(Fig. U: 

St \^ UmVnWl,^{Yi)Wna{Yi)(J-ll^„-ii_„ + h.C. 

(16) 

and 

'^M, 5^U,n\v„^\-^i,„-il,_„ + h.C. (17) 

where UmVn is the coherent factor due to the RVB par- 
ing, with m and n denoting the degenerate lowest energy 
states shown in the inset of Fig. |4] with the degenerate 

Em = Eg- 

It is straightforward to obtain 

X..(r, r) = -i3(-l)"-e"4e-^«- (18) 

and 

Xers(t) ~ ~5De-''^^ (19) 

with D = is the spectral weight whose doping 

dependence is shown in Fig. [6l In Xzz we have used the 

relation^ \j:,^w*mAr)wrnair')\ = . 

Correspondingly the dynamic spin susceptibility is ob- 
tained by 

x'Li<l, ^) = '^-^De-'<^^~^--^'5{uj Eg) (20) 

TT 

and the Aig Raman scattering cross-section 

I^KsH'xx'EKsi^)-5D5{Lo-Eg) (21) 
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FIG. 6: The doping dependence of the spectral weight D for 
the spin-rotons appearing in (|20|l and (|2ip . 



So the triplet spin-roton will appear as a resonancelike 
mode in x"z (Qi ^) uj = Eg, with momentum q peaked 
at the AF wavevector Qaf and a width inversely propor- 
tional to the RVB pairing size Oc which thus determines 
the spin-spin correlation length oc a/VS. Similarly, in 
-fERs('^) a "resonance mode" at Eg will also be exhib- 
ited, which corresponds to the singlet spin-roton excita- 
tion. It should be emphasized that in the neutron and 
Raman scattering measurements only local spinons at the 
same von Neumann lattice are involved and the correc- 
tion from the logarithmic potential Uint in (|13p is always 
negligible. Of course, high-energy spin-roton excitations 
can be also detected by these experiments at w > Eg, 
which will involve spinons at higher energy levels shown 
in the inset of Fig. 01 whose effect^ will not be consid- 
ered in the present work for simplicity. 

At half-filling, the minimal roton energy will be soft- 
ened to zero: i.e., Eg = with Cc ^ oo. As shown 
in Fig. [5] the spectral weight D in ([20)1 remains finite 
at (5 ^ and characterizes the weight of the Neel oder 
as the triplet rotons a,t Eg = are condensed into the 
AF ordering. By contrast, /ers('^) in this limit as 
there is no more charge density fluctuation to couple with 
the incident light in the Raman scattering measurement. 
Furthermore, high-energy triplet spin-roton excitation is 
expected to be reduced to the gapless spin wave^i at 
(5^0 with the spinon spectrum shown in the inset of 
Fig. m becomes a continuous energy spectrum described 
by the Schwinger boson mean-field theory^^. 



C. Tc formula 

We now discuss how thermally excited spin-rotons can 
effectively destroy the phase coherence of the supercon- 
ducting condensation and determine the transition tem- 
perature Tc- 

The long-wavelength superfiuid stiffness ps will be 
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renormalized by spin-roton excitations via the internal 
gauge field A* in the London action ([5|). Such spin- 
rotons shown in Fig. [2] resemble the conventional vortex- 
antivortex pairs in the XY model^, except that the unit 
vorticity of each spinon-vortex is tt instead of 2tt of a 
conventional vortex. A further difference is that the low- 
energy spinon- vortices will distribute on a von Neumann 
lattice with the degeneracy g = 4 as illustrated in Fig. O 
instead of g = 1 on the orginal lattice in the XY model. 
Corresponding to the minimal energy Eg of a spin-roton, 
the fugacity is y = p-Eg/2kBT each spinon effectively 
contributes to a core energy Eg/2. Compared to the 
XY model, such a vortex core energy is much cheaper as 
Eg ^ 5 J at low doping. Thus, the superconducting phase 
transition controlled by spin-rotons, which are governed 
by in (|13p . is expected to be similar to a conventional 
KT transition, but the Tc formula should be quantita- 
tively different due to the peculiar internal structure of 
a spin-roton excitation outlined above. 

In the following, we shall follow a standard textbook 
mathcmatic procedure^ in dealing with a conventional 
KT transition. Define the reduced stiffness K = 



and then the renormalized reduced stiffness Kr, obtained 
by averaging over the spin-roton excitations governed by 
(fT^ . is found by 



Kr = K 



) (22) 



where N is the original total lattice number. The corre- 
lation {qm<lm') can be easily evaluated in terms of 
to lowest order in fugacity y^: 



such that 



=-22/2 [|R^-R^,|/^o] 



K 



R 



K -g^TT^y^K 



dR / R 



Co 



(23) 



(24) 



where the lattice constant of the von Neumann lattice 
provides the short distance cutoff. Again following the 
steps in Ref. [4l|, one arrives at differential renormaliza- 
tion group (RG) equations 



dK- 



dl 



gVy' + 0(y*) 



(25) 
(26) 



with Ku = Kji[K{l), y{l)] remaining as a constant, which 
results in the fixed point at K* = 8/tt and y* — 0. 

Thus, by substituting Kn = limj^tx) ^(0 = ^* i^ito 
p4)) . one gets 



2 fceTe 



(27) 




FIG. 7: [Color online] The coefficient k defined in ([29 
calculated at some typical values of the parameter: t^/J 
and 3, and is weakly doping dependent around k = 6. 



Ps with ~ e ^a/f^BT^^ which can be further rewritten as 



Vc 



8 kBTr. 



(28) 



in which titt with n = 1 denotes the unit vorticity of the 

vortex. (For the sake of comparison, we have introduced 

n in (|28|) such that the case of n = 2 is also allowed 

which corresponds to the conventional 27r vortex in the 

XY model.) Equation ([28|) indicates that the rigidity of 

the superconducting state can only sustain the amount 

9 1 ^ 

of vortex-antivortex pairs with < p- . Using n — 1 
and 3 = 4, one finally finds 



Eg 

ksTc 



21n- 



4^271 



1 - 



which at fc^Tc ^ 7rps/8 results in 

K ~ 21n4\/27r = 5.76 



(29) 



(30) 



Generally, k can be determined self-consistently accord- 
ing to (23) with using nih = l/{2tha'^) and Eg{S) pre- 
sented in Fig. m The result is shown in Fig. [7] as 
a function of doping concentration S at th = 2J and 
th = 3 J, respectively. Fig. [7| indicates that the value of 
K is roughly a universal value at 6 which is weakly de- 
pendent on the choice of as well as the doping concen- 
tration. So we obtain the Tc formula which is in ex- 
cellent agreement with the high-Tc cupratcs as shown by 
the straightline in Fig. [TJ It is noted that y^ = e"" ^ 1 
is consistent with the small fugacity condition used in the 
above derivation of the RG equations. Finally, we com- 
ment that in a previous more complicated approach^, 
Tc was calculated without properly considering both the 
singlet and triplet spin-roton excitations, which resulted 
in somewhat higher and non-universal value of Tc. 
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III. DISCUSSION 

111 this work, we have proposed a consistent under- 
standing of some intriguing experimental facts concern- 
ing high- Tc superconductivity in the cuprates. The key 
concept is the presence of a new type of elementary exci- 
tations in the superconducting state, i.e., spin-rotons, in 
addition to conventional nodal quasiparticle excitations. 
Such novel modes arc composed of supercurrent vortex- 
antivortcx pairs (rotons) locking with free spins at the 
two poles, which form degenerate spin singlet and triplet 
spin states. We have found that they are indeed mea- 
surable by ERS in the Aig channel and INS at the AF 
wavevector as resonancelike modes, which are consistent 
with the experimental observations. In particular, we 
have shown that it is this new kind of excitation that 
determines superconducting phase coherence transition 
with Tc oc -Eg in (U, in excellent agreement with the 
cuprate superconductors. It should be noted that the 
Aig peak has also been probed in the resonant electronic 
Raman scattering experiment^. Similarly, the resonant 
inelastic X-ray scattering (RIXS)^ should be also able 
to detect such a singlet spin-roton excitation if a higher 
resolution (< 40emV) can be achieved. 

So the "resonance energy" Eg as the characteristic en- 
ergy scale of these spin-roton excitations will play an im- 
portant role in the superconducting phase, in contrast to 
the BCS theory in which quasiparticle excitations domi- 
nate. To leading order approximation. Eg vs. doping in 
Fig. m will decide the phase diagram of superconductiv- 
ity. Here Eg (thus Tc) vanishes at overdoping because the 

underlying RVB pairing A* = E^iC^J-ae*"^''^) = 0^^, 
while Eg vanishes at 5 = where the spin-rotons ex- 
perience Bose condensation to form an AF Neel order 
at T = 0. The phase above Tc will be full of free 
spinon-vortices known as the spontaneous vortex phase 
or the lower pseudogap phase^i^^, which may explain the 
Nernst regime discovered^ in the cuprates. 

However, if Eg vanishes at a finite but small doping 
concentration, then the AF order may persist over in a 
finite regime where Tc = 0. As a matter of fact, if the 
non-uniform charge distribution is allowed. Eg as the so- 
lution of PH]) can indeed be softened to zero at some 
small finite doping. A case considering some Z2 topolog- 
ical excitation at low doping does lead to the result that 
Eg vanishes as — Xc at a critical doping Xc — 0.043^. 
Below Xc, either an AF spin glass state or charge stripe 
phases has been show n^^'^^ to be competitive before the 
system becomes a commensurate AF ordered state near 
the half-filling. 

Furthermore, there is no physical reason to protect the 
degeneracy of singlet and triplet spin-roton excitations as 
Eg — > 0. In other words, the residual interaction may de- 
cide which mode will be softened more quickly to result 
in a competing c harg e or spin order at low doping, as con- 
jectured in Refs. [ll| andl48l. The bottomline here is that 
the spin-roton excitations are expected to be essential 
in describing the quantum phase transition of supercon- 



ductivity to other low-doping phases at T = 0. Detail 
investigation along this line is beyond our current scope 
and will be discussed elsewhere. 

Acknowledgments 

We acknowledge helpful discussions with P. W. An- 
derson, D. H. Lee, P. A. Lee, T. Li, N. P. Ong, T. 
Senthil, and X. G. Wen. This work is supported by 
the NSFC grant nos. 10688401, 10834003 and the Na- 
tional Program for Basic Research of MOST grant no. 
2009CB929402. 

APPENDIX A: DIAGONALIZATIONS OF (fTOj) 

The spinon Hamiltonian Hs (fTO|) can be easily 
diagonalized^ under a uniform distribution of the holon 
condensate ph = 6a~^. To be self-contained, in the fol- 
lowing we briefly outline the main results. 

By using the Bogoliubov transformation 

bai'r) = ^{Umlnia " VmlL-a)'^ma{'r) (Al) 
m 

we obtain the spinon Hamiltonian Hg as follows 

Hs^Yl ^rnlLalnia + COUSt. (A2) 

ma 

where 

= .™ = sgn(U)^/5^(A3) 

and 

Em = v/A2 - e2„ (A4) 

Here the quantum number m denotes an eigen-state 
Wraaiji) with the eigen- value 

tmWmaiy^i) ^ -Js ^ C*'^'^?'- W^,^ (r^ ) (A5) 
J=NN(i) 

The spinon excitation spectrum (|A4p exhibits an uneven 
Landau-lcvel-like form as shown in the inset of Fig. |4l 
To obtain this spectrum, we have used a self-consistent 
condition for Js = J(l — 4(5)A''/3^ and the chemical 
potential A in E^ by enforcing bl^ha = (1 - S)N. 

Focusing on the lowest energy level Eg = Eg/ 2, the 
corresponding wave package as the solution of (jASp can 
be express as 

Orn{Vi) = |u7„,^(r,;)|^ ~ cxp j-^jr^ - R„i ^ | 

(A6) 

with the degenerate states labeled^^ by the site R™ in a 
von Neumann lattice shown in Fig. [51 Note that for each 



9 



Rm, there are four degenrate states (g=4) corresponding 
'WmT(r»)> w^ii^i), {-^y'-Wjn^iri) and (-l)'"'^™^ (r^) due 
to the time reversal and bipartite lattice symmetry^. 
Finally one can express the spin operator = 

\ Y.a '^^Icr^i'y in terms of 7,^,^ and 

mncT 
mna 

(A7) 

Here we discard the 77 terms because they have van- 
ishing contribution at the low temperature. And 
for the Raman tensor in the Aig channel^, ta^^ = 
c\^CjrT + h.c, one can use the phase string 
rcprcsentation^fli^ for the electron operator in the t-J 
model and the holon condensation condition to obtain 

TM. = -\Y: /../.le-(^?.+^y6L&,.e-^Uh.c. 

<ij>a 

- E C^.^e^'^^U h.c. 

<ij>a 

cx -(5^u™|w,„|7/„^7^,_^ + h.c. (A8) 

ma 

APPENDIX B: DERIVATION OF Uint IN ([T2II 

According to the discussion in Sec. Ill Al an excited 
spinon will always induce a 7r-vortex as shown in Fig. 



[3l The vortex core will be determined by the spinon 
wavepacket Om{Yi) in (|A6p with a core energy = 
Eg/2. Introduce m = V x A with A = Vcj) + A" to 
describe the winding number for the spinon vortices: 

m(r) = £7r^^O„(ri)(5(r - rj)g„ (Bl) 

where qm{= 0, ±1) denotes the vorticity of a spinon- 
vortex {qm = means no spinon excitation at state m). 
Then by integrating over r in (O in the absence of the 
external electromagnetic field, one can determine an ef- 
fective interaction between the spinon-vortices 

1 f (fa - 
Uint = ^Ps J j^A{q) ■ A(-q) 

1 f d?q m(q) • m(— q) 

= Q^-^Ps^^L - -ps ^ qmqm'Imm' (B2) 

The first term in Uint leads to vortex neutrality Q — 
J2m 1m ~ 0, and in the second term 

imm' = ^Om(rj)0„'(rj)ln|ri - Tjl 




xln|(r^-r;) + (R„-R„OI 
~ ln|(R„-R„OI (B3) 

which leads to 
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